Cosmological tests of general relativity with future tomographic surveys 
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Future weak lensing surveys will map the evolution of matter perturbations and gravitational 
potentials, yielding a new test of general relativity on cosmic scales. They will probe the relations 
between matter overdensities, local curvature, and the Newtonian potential. These relations can be 
modified in alternative gravity theories or by the effects of massive neutrinos or exotic dark energy 
fluids. We introduce two functions of time and scale which account for any such modifications in 
the linear regime. We use a principal component analysis to find the eigenmodes of these functions 
that data will constrain. The number of constrained modes gives a model-independent forecast of 
how many parameters describing deviations from general relativity could be constrained, along with 
w(z). The modes' scale and time dependence tell us which theoretical models will be better tested. 

PACS numbers: 98.62. Sb, 04.80.Cc, 95.80.+p, 98.80.-k 



The observed acceleration of cosmic expansion poses 
a puzzle for modern cosmology. It may be evidence for 
dark energy (DE), a component with a negative equa- 
tion of state, w, that makes it gravitationally repulsive. 
It also warrants studying extensions of general relativity 
(GR) with extra degrees of freedom that can mimic the 
effects of DE. Modifications to GR are well constrained in 
dense regions like our solar system [1]. On larger scales, 
however, GR is less well-tested. Several modifications to 
GR, capable of producing cosmic acceleration have been 
proposed [2]. With the right parameter values, they can 
match the expansion history of a universe made of cold 
dark matter (CDM) and a cosmological constant A - the 
observationally favored ACDM model [3]. However, their 
predictions for the growth of structure can differ since the 
equations for the evolution of perturbations are modified. 
Future tomographic weak lensing surveys, like the Dark 
Energy Survey (DES)[4] and Large Synoptic Survey Tele- 
scope (LSST) [5], will measure lensing shear and galaxy 
counts in many redshift slices (hence the term tomogra- 
phy), thus mapping the evolution of perturbations, and 
offering a new test of GR on cosmological scales [2]. In 
this work, we use a two-dimensional Principal Compo- 
nent Analysis (PCA) to forecast the constraints on mod- 
ified growth (MG) - and thus our understanding of grav- 
ity - coming from these surveys. Unlike previous MG 
forecasts, ours is model-independent, and lets us deter- 
mine how many parameters describing MG could be con- 
strained, along with the regions in parameter space where 
we expect the most sensitivity to MG. 

We consider linear scalar perturbations to the 
flat Friedmann-Robertson- Walker metric in Newtonian 
gauge 

ds 2 = -a 2 (?y)[(l + 2V(x,r]))dri 2 - (1 - 2&{x, r)))dx 2 }, 

where rj is the conformal time and a{j]) the scale factor. 
In Fourier space, one can write [6, 7] 

$/* = 7 (fc,a), fc 2 * = -fj,(k, a)4:7rGa 2 pA , (1) 



where A is the comoving matter density perturbation. 
The function 7 describes anisotropic stresses, while \x de- 
scribes a time- and scale-dependent rescaling of Newton's 
constant G, as well as the effects of DE clustering (a fea- 
ture of many exotic DE models) or massive neutrinos. In 
ACDM, the anisotropic stress due to radiation is negligi- 
ble during matter domination, thus 7 = 1 = fi. In this let- 
ter, we determine how well the unknown functions "f(k, a) 
and fj,(k, a) can both be constrained by future data. We 
also address how well we can detect any departure from 
7 = 1 = /i, without distinguishing between them. 

We consider the two-point correlations (both auto- and 
cross-) between galaxy counts (GC), weak lensing shear 
(WL), and cosmic microwave background (CMB) tem- 
perature anisotropy, plus the CMB E-mode polarization 
and its correlation with the CMB temperature. Detailed 
descriptions of our assumptions for each measurement 
are found in [8] . GC probe the distribution and growth of 
matter overdensities, thus giving an estimate of the New- 
tonian potential $, up to a bias factor. WL is sourced 
by the sum of the potentials + $). GC are also af- 
fected by "magnification bias" [9], where WL conver- 
gence magnifies some faint (thus otherwise undetected) 
galaxies, adding mild dependence on (*F + <I>) to the GC. 
CMB data probe the Integrated Sachs- Wolfe effect (ISW) 
which depends on d(^ + ^)/dr]. Thus, measuring GC and 
WL over multiple redshift bins, along with CMB data, 
yields information about the relation between 'F and 
and their response to matter density fluctuations. Fur- 
thermore, supernovae (SN) redshift-luminosity measure- 
ments and CMB constrain the expansion history. For our 
forecasts, we assume the following probes: Planck [10] for 
CMB, DES and later LSST for GC and WL, and a Su- 
pernova/Acceleration Probe -like the Joint Dark Energy 
Mission (JDEM) [11] for SN. To compare to current data, 
we use large scale structure (LSS) data compiled in [12], 
and the latest SN and CMB data from the 'Constitution' 
SN sample [13] and WMAP5 [3]. For current data, we 
omit WL, but include the GC-CMB cross-correlations. 

In full generality, we treat j(k,a) and /i(k,a) as un- 
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known functions and perform a PCA [14, 15] to deter- 
mine how many of their d.o.f. can be constrained. We 
use redshift z = a -1 — lasa time variable, and pixclizc 
the late-time and large-scale universe (z £ [0, 30], fc £ 
[1CT 5 ,0.2] hMpc" 1 ) into M + 1 z-bins and N fc-bins, 
with each of the (M + 1) X N pixels having independent 
values of fj,ij and jij. We consider w(z) as another un- 
known function, allowing each of the M + l z-bins to have 
an independent value of Wi . Since the surveys we consider 
will not probe z > 3 in detail, we use M bins linear in z 
for z G [0,3] and a single bin for z G [3,30]. We choose 
M = N = 20 and have checked that this pixelization is 
fine enough to ensure the convergence of the results. We 
use logarithmic fc-bins on supcrhorizon scales and linear 
fc-bins on subhorizon scales, to optimize computational 
efficiency. As in [8], we only consider information from 
scales well-described by linear perturbation theory, which 
is only a fraction of the (fc, z)-volume probed by future 
surveys. Since the evolution equations [8] contain time- 
derivatives of jj,(k,z), 7(fc,z) and w(z), we follow [16] 
and use hyperbolic tangent functions to represent steps 
in these functions in the z-direction, while steps in the 
fc-direction are left as step functions. 

Our parameters are the 7^ , [Mj, and w%, along with the 
usual cosmology parameters: energy density of baryons 
Qbh 2 and CDM fl c h 2 , Hubble constant h, optical depth 
r, spectral index n s and amplitude A s . We include one 
bias parameter per GC z-bin, and the intrinsic SN mag- 
nitude. Thus we have (M + 1)(2N + 1) + 17 = 878 pa- 
rameters in total. For a given set of parameter values, 
we use MGCAMB [8, 17], (a modification of CAMB [18] 
developed by us to study modified growth), to compute 
angular spectra for our obscrvablcs. We generate numeri- 
cal derivatives of observables with respect to parameters, 
and use the specifications for the experiments to com- 
pute the Fisher information matrix, which defines the 
sensitivity of the experiments to these parameters (see [8] 
for computational details). Our fiducial values are in all 
cases ACDM: 7^ = /ijj = — Wi = 1 V and the fiducial 
values of the other parameters are those of WMAP5 [3]. 

Let us first study the expected errors on /i(fc, z). The 
error on any /z^ is large, and the pixels have highly corre- 
lated errors. PCA finds the linear combinations of pixels 
with uncorrelated errors. We take only the fiij block of 
the covariance matrix, thus marginalizing over all other 
parameters, including the Wi and 7^. We invert this block 
to obtain the Fisher matrix for our fi values, Fr^), and 
diagonalize i 7 ^) by writing .FL\ = VF T AVF. The rows 
of matrix W are the eigenvectors, or the Principal Com- 
ponents (PC's) [14], while the diagonal elements of A 
are the eigenvalues A m . Each eigenvector, e M (fc,z), is a 
linear combination of the original pixels /iy, forming a 
surface in (fc, z) space. The eigenvectors are orthogonal. 
We normalize them to unity, rescaling the eigenvalues ac- 
cordingly. Then, {e M (fc, z)} forms an orthonormal basis in 
which we can expand /i as //(fc, z) — 1 = Ylm a m e m (k, z), 
where a m are the new uncorrelated parameters with 
variances given by the A m : A m = [cr 2 (o! m )] _1 . We ex- 



pect, from existing data, that variations in fi larger than 
O(l) are unlikely. We enforce this by applying a prior 
A TO > 1 to the matrix F/^y This procedure, analogous 
to the treatment of w(z) in [19], does not affect the well- 
measured modes, but gives a reference point with respect 
to which we define poorly constrained modes. The worst- 
measured modes have variances approaching the prior, 
while those with smaller variances are the well-measured 
ones. Since we compute the full covariance matrix, then 
marginalize over all but the parameter(s) of interest, our 
procedure yields the results that we would get for fi if 
we simultaneously measured w, 7, and /i. This analysis 
can be repeated for 7 or w. Given the PCA results, one 
can convert uncertainties in the expansion parameters 
a. m into uncertainties in any other parameterization of 
fi(-f) without recalculating the Fisher matrices [15, 16]: 
one projects the PC Fisher matrix onto a new basis. 

Measurements probe combinations of $ and ^P, so the 
effects of 7, which affects only are mixed with those of 
H, which affects both potentials. This yields degeneracy 
between fx and 7. By varying both, then marginalizing 
over one, we lose information common to both functions. 
This is necessary when separately constraining fj. and 7. 
While this degeneracy impedes their ability to do so, DES 
and LSST will yield non-trivial constraints on fj, and 7 
with mutual marginahzation (marginalizing over fj, when 
measuring 7, and vice versa). The uncertainties associ- 
ated with PC's of n and 7 are shown in the top two panels 
of Fig. 1 for LSST, DES, and current data. To quan- 
tify the sensitivity of the surveys to MG, we introduce 
three thresholds: well-constrained (Ti, a(a m ) < 0.01), 
constrained (T2, 0.01 < <j{a m ) ^ 0.1), and informative 
(T 3 , 0.1 < a(a m ) < 0.5). From Fig. 1, we see that DES 
could constrain two fi parameters and no 7 parameters. 
LSST could constrain many modes, as it will have a supe- 
rior sky coverage and resolution, wider z-span, and more 
precise photometric redshift measurements. Current data 
effectively cannot constrain cither fi or 7. The constraints 
on /j, are generally stronger than those on 7: [i affects GC, 
WL, and CMB, while 7 primarily affects WL and CMB 
(GC is only affected by 7 via magnification bias). 

Fig. 2 shows selected eigenmodes of // and 7 for LSST 
and DES. The i th fi(-f) mode represents the i th best- 
constrained independent fj,(k, z)( , y(k, z)) surface. Models 
that predict /1 and 7 similar to our "best" modes, then, 
will be better constrained. We observe no degeneracy in 
the fc and z dependences of the modes. This is counter- 
intuitive, since changing fi at some point (fc, z) should 
have the same impact on the observables change 
at a larger scale but later time. However, since we al- 
low for simultaneous variation of 7, the change in fi is 
more efficiently off-set by adjusting 7, eliminating the k- 
z degeneracy. There is a clear pattern to the modes; the 
best constrained modes have no z nodes, but apparent fc 
nodes, and, approximately, the m th mode has m k nodes. 
For LSST, at roughly the 10 th mode, the first z node ap- 
pears, followed by another period of scale-dependent pat- 
terns. The alternating (fc, z) patterns repeat until mixed 
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FIG. 1: Uncertainties in the eigenmodes for current data, and 
for future data sets including LSST(DES). The two upper 
panels show modes of fi and 7 with mutual marginalization. 
The lower panel shows uncertainties in the combined modes. 
The purple and blue solid lines and the shaded region denote 
the thresholds T\,T^ and T3, respectively. The filled symbols 
denote the redshift-dependent modes. In the lower panel, the 
dashed lines show the uncertainties for fi with 7 fixed. 



N 



II 


- | <- 






- - 


>% - 


1 ' 1 ' 1 I. 




1 ' 1 ' 1 ; 

Ml 1 : 


1 ' 1 ' 1 .1 

: M v 


1 ' 1 ' 1 l 


1 ' 1 ' 1 j 


** 

0.05 0.10 0.15 


m 

0.05 0.10 0.15 


0.05 0.10 0.15 


- <- 

0.05 0.10 0.15 



k h/Mpc 

FIG. 2: Eigensurfaces for and 7, with mutual marginaliza- 
tion, for LSST(DES) along with Planck and JDEM. 



(k, z) modes appear, after which there are no clear pat- 
terns. The best modes are mainly functions of k and not 
z. This is partly because the total observable volume in 
the radial (z) direction is limited by the dimming of dis- 
tant objects and, ultimately, the fact that structures only 
exist at relatively low z. Also, it is related to us consid- 
ering only linear perturbations in our analysis, since at 
small z the observable volume is too small to fit the small 
fc-modes that are still in the linear regime. Hence, there is 
more volume available for studying the spatial distribu- 
tion of structure than the radial distribution. The num- 
ber of nodes in the z and k directions tell us, respectively, 
the number of z- and fc-dependent parameters that sur- 



veys could constrain. For example, for LSST, there are 
three clear z-dependent patterns of \i and 7 whose eigen- 
values fall within T2 and T3 ranges. For DES, the only 
constrained modes are those of /x, and exhibit only one 
type of z-dependence, with one node. 

The LSST modes for [i and 7 have similar shapes ex- 
cept that the fj, modes have a deeper z span than the 7 
modes. This is due to an accumulation effect on fi. On 
subhorizon scales, the density contrast, (related to \& via 
Eq. (1)), evolves via A + HA = 4-KG(j>pa 2 A. Perturbing 
fi at one pixel, e.g. enhancing it at some k and z = 3, 
enhances A (and thus \&) for that fc-mode V z < 3, since 
the growth factor at all later times is enhanced. On the 
other hand, changes in 7 at high 2, which primarily affect 
WL through changes in $ (Eq. (1)), do not affect WL at 
low z. Hence, the z-sensitivity range of 7 is primarily de- 
termined by the redshift range of the WL kernel. The 
z-dependence of 7 also affects the ISW contribution to 
the CMB at small k and low z, but its contribution to 
the Fisher matrix is small due to a large cosmic variance. 

In most models of modified gravity, fj, and/or 7 evolve 
in a time- and scale-dependent way [2]. For example, in 
scalar-tensor theories, they undergo a step-like transition 
at the Compton scale of the scalar field. The peaks in 
the eigensurfaces indicate the "sweet spots" in the (k, z) 
space where such a transition scale can be detected by 
the survey, while the frequencies of the modes tell us 
how well a transition can be resolved. For LSST, the 
best-measured modes peak in the region 0.04 < k < 
0.16h _1 Mpc and 0.5 < z < 2, indicating a sensitiv- 
ity to a Compton scale today of 50 < A° < 1500 Mpc, 
where we have allowed for a range of possible time evo- 
lutions of the mass scale [8]. Massive neutrinos also intro- 
duce a transition in /j, due to free-streaming. From the 
expression for the free streaming length in terms of z 
and the neutrino mass m„ (see e. g. [22]), we find that 
the transition scale is within the LSST sensitivity window 
for 0.1 < m u < 0.7 cV. Smaller masses induce an over- 
all suppression of growth with no scale-dependent signa- 
tures. While observable to some extent, this suppression 
is largely degenerate with w(z), especially if one allows 
for an arbitrary evolution of w(z), as we have done. 

In addition to constraining \i or 7 individually, a less 
ambitious yet equally interesting question is how sensitive 
data is to any departure from standard growth. Namely, 
one may ask if either function deviates from unity, with- 
out specifying which. For this purpose, we want to save 
the information common to both functions, which we pre- 
viously lost by mutual marginalization. Hence, we con- 
sider the combined principal components of \i and 7. We 
follow the same procedure as before, except now we di- 
agonalize the block of the Fisher matrix containing (i 
and 7 pixels. The eigenvalues of these combined PC's are 
shown in the lower panel of Fig. 1. Even today's data can 
provide around 15 "constrained" modes. This is unsur- 
prising since the LSS power spectrum, P(k), is known to 
better than 10% precision over an order of magnitude in 
k. Changing fi at any k directly affects P(k), which means 
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that large variations in /i arc disallowed in the range 
where P(k) is well-measured. After marginalizing over 7, 
the direct impact on P(k) is lost, since one can now offset 
changes in \i by adjusting 7. The dashed lines in Fig. 1 
represent the eigenvalues of fi without marginalizing over 
7. They are comparable to the eigenvalues of combined 
PC's, supporting the above notion. From Fig. 1, we see 
that DES will provide around 20, while LSST will provide 
up to 100, "well-constrained" combined modes. 

So far, our analysis has neglected systematic errors, 
which are model-dependent and very hard to predict. 
While we address systematics thoroughly in an upcoming 
publication, we report here the outcome of a preliminary 
analysis based on the assumptions in [20, 21]. We repeat 
our PCA with extra parameters describing likely sources 
of systematic error: shifts in the centroids of the z-bins, 
distortions of the z-bin distribution functions, and addi- 
tive and multiplicative errors on the WL signal due to 
point-spread- function contributions, as in [20] . This adds 
58(80) parameters to our analysis for DES(LSST). We as- 
sume no "catastrophic" photo-z mis-estimation, and ap- 
ply a conservative set of priors [21] to these parameters 
and marginalize over them. We find that the systemat- 
ics result in a noticeable, but not a dramatic, dilution of 
constraints on MG from DES. This is because photo-z 
errors would most immediately affect the z-dependence 
of MG, to which DES was only weakly sensitive even 
without the systematics. As discussed above, constraints 
from DES will be primarily on the scale-dependence of /x 
and 7, and that information is mostly preserved. The im- 
pact of the systematics on LSST forecasts is more signif- 
icant, because LSST has a higher potential for resolving 
z-dependent features. There too we find that inclusion of 
systematic errors preserves most of the scale-dependent 
information but can reduce our ability to measure eigen- 
modes of /1 with z-dependent features, underscoring the 
need to study and control systematic errors in lensing 
surveys. Even after accounting for systematics, LSST and 
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DES are powerful probes of MG. 

The Dark Energy Task Force [23] recently analyzed 
constraints on DE from future surveys, without consid- 
ering MG, and found that w{z = 0) and (dw/dz)\ z =o 
could both be constrained. A time-varying w(z) alters 
the growth dynamics in a scale-independent way, so the 
scale-dependence of fj, and 7 cannot be duplicated by a 
choice of w. Furthermore, since we consider linear scales, 
the dominant portion of the information on MG comes 
from higher z (z > 0.5), at which DE effects are not as 
important. Thus, in addition to measuring w(z), future 
surveys will tightly constrain scale-dependent departures 
from ACDM, and those occurring at high redshifts. Note 
that including the non-linear growth data from lower z 
requires a model dependent treatment of MG, in which 
case w(z) and MG would be related by the same theory. 

To recap, we have forecasted the constraints on mod- 
ifications to GR from future data sets, in comparison 
with existing data. We find that combined data from 
Planck, JDEM and LSST(DES) can tightly constrain 
around 100(20) parameters of MG (corresponding to the 
100(20) eigenvalues in region Tl in the combined fi and 7 
analysis), if the systematics are negligible. Current data 
can constrain only one parameter to this level. We have 
further identified the regions in parameter space to which 
future datasets are most sensitive. In general, our tech- 
nique can be used in survey design to move the sweet 
spots to the most interesting parts of parameter space. 
Our results are obtained using only linear-scale data, be- 
ing a conservative "proof of concept" that upcoming sur- 
veys can rigorously test GR over cosmic distances. 
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